Critical Rotational Frequency for Superfiuid Fermionic Gases across Feshbach 

Resonance 
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We present a method to determine the critical rotational frequencies for superfluidity of both 
uniform and trapped Fermi gases across wide Feshbach resonance. It is found that as one approaches 
the resonance from the BCS side, beyond a critical scattering length, pairing is so robust that 
superfluidity cannot be destroyed by rotation. Moreover, the critical frequency has a sequence of 
jumps revealing the appearance of Landau levels, which are particularly prominent for systems 
up to a few thousand particles. For rotational frequency below an "ultimate" critical frequency, 
defined to be the lowest frequency at which the center of the cloud goes normal, a trapped gas has 
a superfiuid core surrounded by a normal gas as seen in recent experiments [C. H. Schunck et.al. 
|cond-mat/0607298l . 



The recent discovery of vortex lattice in a rotating 
Fermi gas of 6 Li near Feshbach resonance by the MIT 
group [HQ provides first direct evidence of phase coher- 
ence of this dilute system. This is analogous to the vor- 
tex lattice in type-II superconductors in the presence of 
magnetic field. Like superconductors, which turn normal 
in sufficiently large magnetic fields, fermion superfluids 
will also turn normal for sufficiently large rotational fre- 
quency. This question of critical rotational frequency be- 
comes even more interesting at resonance since fermion 
superfluidity is strongest at resonance and the properties 
of the system is universal in the absence of rotation^. 
As rotation increases, more angular momentum is de- 
posited into the system, driving it toward the quantum 
Hall regime. How the fermion system copes with the 
strong pairing interaction and large angular momentum 
is an intriguing and fundamental question. 

For harmonic traps, the external rotation frequency Q 
must be less than the radial frequency of the harmonic 
trap to in order to confine the system. The angular mo- 
mentum of the system increases as f2 increases, and be- 
comes a maximum when Q/lu = 1. Since critical rotation 
frequency reflects the strength of superfluidity, it scales 
with fermion density and is therefore a decreasing func- 
tion of radial distance, il c = O c (r) . Thus, as f2 increases, 
the surface of the cloud will first turn normal, while the 
center region remains superfiuid with an array of vor- 
tices. At a particular frequency f2*, the center region of 
the superfiuid shrinks to zero and the whole cloud turns 
normal. It is this "utimate " critical frequency at very low 
temperatures that is the subject of our discussion. 

In this letter, we shall derive the phase diagram of the 
rotating fermions at temperatures kgT << huj and dis- 
cuss its unique features. The result is shown in fig[T] 
We find that: (1) As one approaches the resonance from 
the BCS side, Q* increases and has a sequence of steps 
which reflect the appearance of Landau levels. The size 
of these steps shrinks as particle number increases, but 
becomes sizable for systems with a few thousand parti- 
cles. (2) fl*/uj remains unity until the scattering length is 



beyond a critical value a* that depends on particle num- 
ber, (see figH]). This means the strong interaction near 
resonance is able to protect pairing even in the limit of 
large angular momentum. While this may not be obvi- 
ous, it can understood by considering the opposite. If 
fl* < lu at resonance, the system would be normal when 
n approaches to and all fermions would be in the lowest 
Landau level. As one makes the system more and more 
dilute, one eventually reaches a two-body problem. How- 
ever, explicit solution of two-body problem shows that at 
resonance, the fermions do not stay at the lowest Lan- 
dau level, as higher Landau levels enable the fermions to 
come closer to take advantage of the strong interaction. 
To eliminate this contradiction, the ratio tt*/u) can only 
be less than 1 on the BCS side of the resonance. 

(A). How f2* is determined: The Hamiltonian of 
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FIG. 1: fij is the "ultimate" critical frequency defined as the 
critical frequency at the center of the trap. a s is scattering 
length. This diagram is for a system of 10 3 fermions. The 
aspect ratio of the trap is uj/uj z — 2.5, where uj and lo z are 
radial and axial frequencies, d = ^Jh/muj. The temperature 
range is k B T/huj = 10~ 3 . (SF) and (N) denote superfiuid and 
normal phase respectively. The insert is a zoom-in plot of the 
strongly interacting regime. 
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a Fermi gas in a rotating harmonic trap with chemical 
potential p is K = H — VlL z — pN = 7Y(fi) — J p(r)n(r), 
where W(fi) = ^/li(fi) + YU>j U v is tne many-body 
Hamiltonian of a uniform system with single particle 
Hamiltonian /i(fi) = (p — mfiz x r) 2 /(2m), and t/y is the 
interaction between particles. The quantity p(r) is de- 
fined as p(r) = fi— V(r), where V(r) = m(u} 2 — Q 2 )r 2 ± /2+ 
mui 2 z 2 /2, r± = (x, y). Our strategy is to first find fi c for 
the uniform system /C (O, p) = 7i(fi) — pN, and then use 
local density approximation to study the trapped case. 

As we shall show later, the condition for superfluidity 
for a uniform system with chemical potential p is that 
the pairing susceptibility W (a function of p/(Ml) and 
defined later) to be sufficiently large, satisfying 

-a/a a <W(n/m), a = y/h/2mfl = a(fi). (1) 

The critical frequency fi c is determined by eq. (fT]) with an 
equal sign. The number density of this uniform system 
is denoted as n, and n = n(p, fi). 

The situation in a trap is trickier, since in that case the 
chemical potential and density at r is /i(r) = p — V(r) 
and n(r) = n(//(r),fi) respectively, p itself is deter- 
mined by the number equation TV = Jn(/x(r), O), which 
gives N = N(p,fl). To find the critical frequency at 
r, one inverts the number relation to get p = p(N, fi), 
and then solves the equation — a(fi)/a s = W([p(N, fi) — 
V(r)]/Ml), which gives fi c (r). In general, when the sys- 
tem has a superfluid core surrounded by a normal cloud, 
the equation N = J n(r) is quite complex for it requires 
the knowledge of n of a uniform gas in both normal and 
superfluid phase. The situation is simplified considerably 
for the ultimate critical frequency fi*,. In that case, the 
whole cloud is normal, only the center is barely super- 
fluid. As a result, n(r) is simply the density n N (p(r)) of 
a normal gas in a rotating trap, where n.jv(/x) is the num- 
ber density of a uniform rotating normal gas, which is 
much easier to calculate; and N = J n N (r) = N N (fi, fi). 
The corresponding expression for p after inverting this 
equation is p = p N (N, fi). The equation for fi* is then 

where d = yj h/ (muj). The solutions of eq.© give figH] . 

(Bl) Solving 7Y(fi) in the Landau gauge: For uni- 
form systems, it is convenient to express the eigenstates 
of /i(fi) in terms of those in Landau gauge: <f> v {x,y, z) = 
e ikz e imQ.xy/h e -ipx Un ^y^ wnere v = ( nj p;fc) labels the 

Landau level n and degeneracy p, and momentum k along 
z; u np = H n (y/a- pa)e-(y/ a -^ 2 / 2 /y/2^^, H n is 
the Hermite polynomial with index n. We have applied 
periodic boundary condition along z and x, taking the 
sample size L z and L x to be 1. We also find it is conve- 
nient to use two-channel model description. The eigen- 
states for close channel molecules have the same form 
as <t> v , except that m is replaced by 2m and a is replaced 



by a/y2. It is also sufficient to write the two-channel 
model in the rotating frame as, 

va vv' 

(3) 

where f£ a creates an open channel fermion with quan- 
tum number v, spin a, and energy = h 2 k 2 /(2m) + 
(2n + I) Ml — p. (£„ is independent of p). creates 
a close channel molecule. Since we shall focus on the 
ground state, and since the condensate occurs at the 
lowest energy, it is sufficient to restrict the Bose field 
b to in the lowest Landau level and to zero momen- 
tum along z|4(, with energy So = 7 + Ml — 2p, where 
7 is the unrenormalized detuning between close chan- 
nel molecule and open channel fcrmions. It is impor- 
tant to note that restricting the condensate to the low- 
est energy state imposes no restrictions on the quantum 
numbers of the fcrmions. Two fermions in different Lan- 
dau levels can convert into a molecule in the lowest en- 
ergy state, with the coupling a and an overlap integral 
Q vv , = J drdR0*(R + r/2)0*,(R-r/2)L>(R,r), where 
R and r are the center-of-mass and relative coordinate 
of the Fermi pair, and D(R, r) is the wavefunction of the 
tightly bound close-channel molecule, which can be taken 
as $(R)<5(r). The field ^(R) has zero momentum along 
z and is a general superposition of in the lowest Landau 
level, i.e. $(R) = J2 q c «3^(o,g,o)j where c 9 's are arbitrary 
except for the normalization X! 9 l c q| 2 = An explicit 
calculation in Appendix I gives 

_ (-1)" H N (V2pa)e-^ 2 

where p — p — p' , P = p + p' , and N = n + n' . For later 
use, we also note that 

E \Q»" 1 2 = ^ n > N ) E m 2 = ^ n > *o> ( 5 ) 
p,p p 

1 AH 
^^ = 23/W2^n!(jV- n )! - (6) 

(B2) Renormalzation: To solve the resonance 
model, it is essential to relate the parameters (a, 7) 
to physical properties such as scattering length a s and 
effective range r . This is done by solving the two 
body problem. The procedure is explicitly shown in Ap- 
pendix II. The result is that 7 is formally divergent, and 
must be renormalized to a finite parameter as 7 = 7 
— a 2 m/(4:Tr\/2ah 2 ) J2 N 1/y/N + 1. The finite parameters 
(7, a) must be chosen to produce the low energy scatter- 
ing properties (a s ,r ). The relation are 

_ J_ _ m J_ _ r Q m 2 

a 2 ~ a s h 2 TrV2' a 2 ~ 2ttV2H 4 ' 
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(CI) Critical Frequency for uniform system. For 

the many-body case, within mean-field approximation 
a(b) = A, H(Ci) becomes 



n 



(8) 

Near tt c , A is sufficiently small. The free energy T ex- 
panded in powers of A is 

_ SpA 2 v tanh(^) + tanh(%i) ,2-,^- 

« 2 ^ m+&) lQl "" ] A+tAA 

(9) 

where /3 = l/(fcsT) and the coefficient in front of A 4 
is positive. Minimizing T with respect to A, using the 
renormalization condition for 7 and eq.([7]), and at the 
end using the condition for wide resonance r Q J a — ► 0, we 
obtain the linearized gap equation, eq.([T]), 

a _ ir\/2ah 2 tanh(^) + tanh(^) 2 



fi k B T 



Ml 



(10) 



Solution of Eq. fTU)) gives the critical frequency tt c as a 
function of \i and temperature T. 

Due to the constraint on Q vv t in 1^, the sum Ylw' 
is over (p, P\ n, n'; fc), and + £ u i = h 2 k 2 /m + (27V + 
l)Ml — 2// = E^k — 2^i. If [i does not coincide with 
any Landau level, those Laudau levels below \i will lead 
a vanishing energy denominator for some non-zero value 
of k, which in turn causes the fc-sum to be logarithmic 
divergent if it were not the ian/i-factors in eq. (|10|) which 
eliminate this divergence at non-zero temperatures. On 
the other hand, if \i coincides with a Landua level N, the 
corresponding energy denominator vanishes at k = 0, the 
fc-sum becomes J dk/k 2 and is divergent. This will lead 
to a sequence of spikes in W as a function of /i. 

To express ea. (|10l) in a form convenient for calcula- 
tion at low temperatures, we note that W is made up of 
the following contributions: (A) both (£„ , £„/ ) have the 
same sign, (B) both (£i/,£zv) have opposite signs, corre- 
sponding to fermions on the same side or opposite sides 
of the chemical potential respectively. Contribution (A) 
can be further divided into: (Al) \£ v ,\ > k B T, (A2) 
l^'l < ^bT. In the case of (Al), the tank-terms in 
eq. dTU)) can be replaced by 1. On the other hand, the 
tanh-terms make the summand of (A2) and (B) finite, 
and therefore become insignificant as T — > 0. One can 
therefore approximate W for k B T « Ml as 



k R T 



m mi 



iry/2ah 2 



E 



1 



where e(x) = ±1 if x<0. Note that both and £ u i are 
independent of p and p' 7 we can therefore sum over the 
degenerate index P and p in eq. lfTTj) as in the two atom 
case and obtain a result independent of {cp}. This shows 
n c is insensitive to the configurations of vortices at their 
onset in the superfluid phase, as in type-II superconduc- 
tor. We then obtain 



E 

N 



4VA+T 



(11) 



ah 2 



dk£(N, k) 



1 



2^ 2 mJ Q [{E Nk -2p) 2 + 8 2 y/ 2 



C(n, k) 



N 

E 

n=0 



£{$v) + £{£.w)\n n , N , 



(12) 

(13) 



where v — (n,p;k) and v' = (N — n,p',—k). While 
many more simplification of can be made, which we shall 
not present for length reasons, eq. (fT2"|) and (fT3"|) are in 
a form convenient for numerical calculations. We have 
calculated W for a system of 10 3 fermions from eq. (fT2"|) 
for k B T<lO~ 3 H£l and have verified numerically that W 
changes insignificantly even when T is lowered by another 
factor of 10. Thus, from now on, we shall not display the 
T dependence of W explicitly. 

To determine the critical frequency f2 c , we rewrite 
eq. (fT0f as -l/(k a s ) = x~ 1 / 2 W(x) = Z(x), where 
h 2 k 2 /(2m) = fx, x = (j,/ (Ml). Fig^A) displays the solu- 
tions of this equation, which are the interactions between 
Z(x) and the horizontal line y = —l/(k a s ). The spikes 
in Z(x) is due to the presence of Landau levels. The 
regions of x where Z > —l/(k a s ) and Z < —l/(k a s ) 
correspond to superfluid and normal phase respectively. 
The largest intersection x* corresponds to the lowest ro- 
tational frequency at which the superfluid turns normal. 
Hence, we identify fl c = (fi/h)/x* . Since Z(x) is positive 
definite, eq. (fTU| only has a solution if a s < 0, i.e. on 
the BCS side of the resonance. Moreover, since Z has a 
minimum (located at x — x a — 2.67), eq. ljlOp has no so- 
lution for k a s > k a a* s = —1/Z(x ), Z(x ) = 0.9716. 
This means approaching the resonance from the BCS 
side, the critical frequency will rise to Q Q = (p/h)/x 
at k Q a* s = —1/Z(x ) = —1.029. For scattering length 
a s > a*, the system will remain supefluid for all f2. 

(C2) "Ultimate" critical frequency O* in a trap. 
The scheme of finding fi* has been discussed in (A) . We 
shall start with the number density of a uniform rotating 
normal system, which is njv(/z, Vl) — 2 ^2 npk Q(fi— (2n + 
l)Ml - h 2 k 2 /(2m)), where Q(x) = 1(0) if x > 0(< 0). 
Summing over k and using local density approximation, 
we have 



n N (fi(r),n) = 



1 



Mr) 
Ml 



-2n- 1, 



(14) 



where n* = [/x(r) / (2Mi) — 1/2] and [x] means the largest 
integer smaller than x. Eq. (fT4")) allows one to find the 
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FIG. 2: (A)Z(x) for a uniform system, x = (j,/{hQ). As one 
approaches resonance from the BCS side, the horizontal line 
comes down. The last intersection takes place at scattering a* 
when the horizontal line touches the minimum of Z(x). (B) 
The function G(Q/lo) plotted for No = 10 3 fermions and the 
aspect radio uj/lo z = 2.5. (C) The function of Y(r/d) plotted 
for n = Whu, u)/u) z — 2.5, fl/cu = 0.8 and z — 0. 



relation between /i and N in the rotating normal state, 
and to construct the function G{N,Q). In FigOjB), wc 
have shown the function G{N, fi] for an anisotropic trap 
(uj/uj z — 2.5) (like that in rcf. l|) with the scaled vari- 
able il/u>. The ultimate critical frequency fi* is given 
by the intersection between G(fi/w) and the horizontal 
line —d/a s with smallest fi. It is from these intersec- 
tions that figH] is constructed. It should be noted that 
fi/w can at most be 1, and G(fi/w) is positive definite 
with a global minimum at fi*/w close to 1. The situa- 
tion is almost identical to (CI): fi* exists only on the 
BCS side of the resonance. It undergoes a series of jumps 
as one approaches the resonance, reflecting the presence 
of Landau levels. It will finally reaches fi* when the 
scattering length reaches a critical value a* such that 
—d/a* = G(fi*/w). For all a s > a*, the system remains 
a superfluid even when fi/w reaches 1. Calculations with 
larger number of atoms show that the phase boundary in 
figQ]will be pushed up to higher rotating frequency and 
the size of Landau step will shrink. However, the basic 
features will still hold. For systems with few thousand 
particles, (see fig.l), the Landau steps are highly visible. 

(C3) Separation of vortex and normal regions: 
Eq. © also provides a convenient means to determine the 
spatial distribution between vortex and normal regions 
as found in the MIT experiment Q , as the superfluid re- 
gion must have a sufficiently large pairing susceptibil- 
ity. Rewriting eq.© in dimensionless form, and in lo- 
cal density approximation, condition for superfluidity at 
r is -d/a s < y/2ST/wW(p-V(r)/(m)) = Y(r/d). In 



fig. 2c, we have plotted Y(r) as a function of r. The region 
of superfluid can be identified easily from the intersection 
between Y and the horizontal line y = — d/a s ^. 

We have thus established the results (1) and (2) 
mentioned in the beginning. Finally we would like to 
point out that although we consider the case where 
&bT/(M7) ~ 10~ 3 , our results should hold as long as 
fesT/(Sfi) < 1. This condition is satisfied in traps with 
frequencies about 300 Hz (~ 10~ 8 K) at the lowest tem- 
perature attainable today (10 _9 K). We would also like to 
point out that although mean field theory is not exact, it 
is known that it provides a good qualitative description 
at very low temperatures. Fluctuations effects at tem- 
peratures higher than Ml will be considered elsewehre. 

Appendix I. The coefficient Q vv ' is given 

by E^/rf 3 RC(R)^(R)<M R), which is 

Q»v' = E q CqUu'KP + P'- g)/V7r 3 / 2 a2(™+™'- 1 )n!n'!, 
denoted as (*), where l vv i = J^°° dtH n (t — t )H n >{t + 
io)e~ t2 e~( t2+t o), to — (p — p')a, and t = y/a. Using 
the well known relations H n (t) = {-l) n e t2 d n (e- t2 ) / dt n 
and e -s 2 +2(t-t ) s _ Yl=o H n{t - t )s n /n\, integra- 
tion by parts and straightforward calculations gives 
hv> = (-l) , VV2 (n+n ' +1) ^n+n'(V2to)e-*«. Substitut- 
ing this result back to eq.(*) gives eq.((4]). 

Appendix II. The two-body eigenstate with en- 
ergy E is = [Btf + ^ A^f^ftjO), where 
(E - j)B - Q UV >A VV >, A vpl {k) = aQ* vu ,B/(E - 

E Nk ), where E Nk = 2(N + + h 2 k 2 /m. 

Hence, we have E - 7 = a 2 Y^ vu > \Q^'\ 2 /( E - E Nk ). 
Eq.©-© imply E„,p, ? IQ-'| 2 = En=o v(n,N) = 
l/(2 3 / 2 7ra 2 ). The equation for bound states is 7 — E = 
a 2 m/{ATT^/2ah 2 ) ~£ N ±/V N + 1 - E/{Ml). The diver- 
gence in the sum can be removed using the renormal- 
ized detuning 7 defined in section (B2), and we have 
7 - E = a 2 m/{4irV2ah 2 )J2 N (l/\/N +1- E / {Ml) - 
1/y/N + 1). To relate (7, a) to {a s ,r ), we take the limit 
fifi -> 0, which yields 7 - E = -a 2 m 3 / 2 % / zr E/(7r v / 2/i 3 ). 
Comparing it to equation derived from scattering theory, 
y/—E = h/ {a s y/m) — roy/mE / {2h), we find eq.©. 
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